Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe.

BACKGROUND: Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry. METHODS AND RESULTS: MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank. CONCLUSIONS: Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

A chromosome-level genome assembly of the yellow-throated marten (Martes flavigula).

The yellow-throated marten (Martes flavigula) is a medium-sized carnivore that is widely distributed across much of Asia and occupies an extensive variety of habitats. We reported a high-quality genome assembly of this organism that was generated using Oxford Nanopore and Hi-C technologies. The final genome sequences contained 215 contigs with a total size of 2,449.15 Mb and a contig N50 length of 68.60 Mb. Using Hi-C analysis, 2,419.20 Mb (98.78%) of the assembled sequences were anchored onto 21 linkage groups. Merqury evaluation suggested that the genome was 94.95% complete with a QV value of 43.75. Additionally, the genome was found to comprise approximately 39.74% repeat sequences, of which long interspersed elements (LINE) that accounted for 26.13% of the entire genome, were the most abundant. Of the 20,464 protein-coding genes, prediction and functional annotation was successfully performed for 20,322 (99.31%) genes. The high-quality, chromosome-level genome of the marten reported in this study will serve as a reference for future studies on genetic diversity, evolution, and conservation biology.

A novel vesivirus (family Caliciviridae) in European badgers (Meles meles) in Hungary, 2020/2021.

In this study, a novel vesivirus (family Caliciviridae) was detected and characterized in faecal and tissue (blood and spleen) specimens collected from three (23.1%) out of 13 European badgers (Meles meles) in Hungary that were tested using RT-PCR and sequencing methods. The complete genome of the vesivirus strain European badger/B40/2021/HUN (OQ161773) is 8,375 nucleotides in length. The ORF1, ORF2, and ORF3 proteins have 81.1%, 70.5%, and 64.2% amino acid sequence identity, respectively, to the corresponding proteins of Asian badger vesivirus, which was first reported in badgers in China in 2022. These results indicate that more than one lineage/species of vesiviruses circulates in mustelid badgers in geographically different regions.

Maps of Constitutive-Heterochromatin Distribution for Four Martes Species (Mustelidae, Carnivora, Mammalia) Show the Formative Role of Macrosatellite Repeats in Interspecific Variation of Chromosome Structure.

Constitutive-heterochromatin placement in the genome affects chromosome structure by occupying centromeric areas and forming large blocks. To investigate the basis for heterochromatin variation in the genome, we chose a group of species with a conserved euchromatin part: the genus Martes [stone marten (M. foina, 2n = 38), sable (M. zibellina, 2n = 38), pine marten (M. martes, 2n = 38), and yellow-throated marten (M. flavigula, 2n = 40)]. We mined the stone marten genome for the most abundant tandem repeats and selected the top 11 macrosatellite repetitive sequences. Fluorescent in situ hybridization revealed distributions of the tandemly repeated sequences (macrosatellites, telomeric repeats, and ribosomal DNA). We next characterized the AT/GC content of constitutive heterochromatin by CDAG (Chromomycin A3-DAPI-after G-banding). The euchromatin conservatism was shown by comparative chromosome painting with stone marten probes in newly built maps of the sable and pine marten. Thus, for the four Martes species, we mapped three different types of tandemly repeated sequences critical for chromosome structure. Most macrosatellites are shared by the four species with individual patterns of amplification. Some macrosatellites are specific to a species, autosomes, or the X chromosome. The variation of core macrosatellites and their prevalence in a genome are responsible for the species-specific variation of the heterochromatic blocks.

Blue Turns to Gray: Paleogenomic Insights into the Evolutionary History and Extinction of the Blue Antelope (Hippotragus leucophaeus).

The blue antelope (Hippotragus leucophaeus) is the only large African mammal species to have become extinct in historical times, yet no nuclear genomic information is available for this species. A recent study showed that many alleged blue antelope museum specimens are either roan (Hippotragus equinus) or sable (Hippotragus niger) antelopes, further reducing the possibilities for obtaining genomic information for this extinct species. While the blue antelope has a rich fossil record from South Africa, climatic conditions in the region are generally unfavorable to the preservation of ancient DNA. Nevertheless, we recovered two blue antelope draft genomes, one at 3.4× mean coverage from a historical specimen (∼200 years old) and one at 2.1× mean coverage from a fossil specimen dating to 9,800-9,300 cal years BP, making it currently the oldest paleogenome from Africa. Phylogenomic analyses show that blue and sable antelope are sister species, confirming previous mitogenomic results, and demonstrate ancient gene flow from roan into blue antelope. We show that blue antelope genomic diversity was much lower than in roan and sable antelope, indicative of a low population size since at least the early Holocene. This supports observations from the fossil record documenting major decreases in the abundance of blue antelope after the Pleistocene-Holocene transition. Finally, the persistence of this species throughout the Holocene despite low population size suggests that colonial-era human impact was likely the decisive factor in the blue antelope's extinction.

Protection status, human disturbance, snow cover and trapping drive density of a declining wolverine population in the Canadian Rocky Mountains.

Protected areas are important in species conservation, but high rates of human-caused mortality outside their borders and increasing popularity for recreation can negatively affect wildlife populations. We quantified wolverine (Gulo gulo) population trends from 2011 to 2020 in > 14,000 km2 protected and non-protected habitat in southwestern Canada. We conducted wolverine and multi-species surveys using non-invasive DNA and remote camera-based methods. We developed Bayesian integrated models combining spatial capture-recapture data of marked and unmarked individuals with occupancy data. Wolverine density and occupancy declined by 39%, with an annual population growth rate of 0.925. Density within protected areas was 3 times higher than outside and declined between 2011 (3.6 wolverines/1000 km2) and 2020 (2.1 wolverines/1000 km2). Wolverine density and detection probability increased with snow cover and decreased near development. Detection probability also decreased with human recreational activity. The annual harvest rate of ≥ 13% was above the maximum sustainable rate. We conclude that humans negatively affected the population through direct mortality, sub-lethal effects and habitat impacts. Our study exemplifies the need to monitor population trends for species at risk-within and between protected areas-as steep declines can occur unnoticed if key conservation concerns are not identified and addressed.

Badger genome grows Tree of Life.

Sequencing the genome of the European badger is the latest step in a project to determine the genomes of all species in the UK and beyond. Georgina Mills reports.

Identification of mutant gene for Black crystal coat and non-allelic gene interactions in Neogale vison.

Sable (Martes zibellina) and American mink (Neogale vison) are valuable species characterized by a variety of coat colour produced on fur farms. Black crystal fur phenotype is Mendelian codominant trait: heterozygous animals (Cr/ +) have white guard hairs scattered predominantly on the spine and the head, while homozygous (Cr/Cr) minks have coats resembling the Himalayan (ch/ch) or white Hedlund (h/h) types. It is one of the most recent of more than 35 currently known phenotypic traits of fur colour in American mink. Black crystal fur phenotype was first described in 1984 in the Russian population of mink, which had undergone selection for domestic defensive response to humans. Here, we performed whole-genome sequencing of American mink with Cr/Cr phenotype. We identified a missense mutation in the gene encoding the α-COP subunit of the COPI complex (COPA). The COPI complex mediates retrograde trafficking from the Golgi system to the endoplasmic reticulum and sorting of transmembrane proteins. We observed an interaction between a newly identified mutation in the COPA gene and a mutation in the microphthalmia-associated transcription factor (MITF), the latter mutation led to the formation of the white Hedlund (h/h) phenotype. Double heterozygotes for these mutations have an entirely white coat and a black-eyed phenotype similar to the phenotype of Cr/Cr or h/h minks. Our data could be useful for tracking economically valuable fur traits in mink breeding programs to contribute to global fur production.

Chromosomal-level reference genome assembly of the North American wolverine (Gulo gulo luscus): a resource for conservation genomics.

We report a chromosomal-level genome assembly of a male North American wolverine (Gulo gulo luscus) from the Kugluktuk region of Nunavut, Canada. The genome was assembled directly from long-reads, comprising: 758 contigs with a contig N50 of 36.6 Mb; contig L50 of 20; base count of 2.39 Gb; and a near complete representation (99.98%) of the BUSCO 5.2.2 set of 9,226 genes. A presumptive chromosomal-level assembly was generated by scaffolding against two chromosomal-level Mustelidae reference genomes, the ermine and the Eurasian river otter, to derive a final scaffold N50 of 144.0 Mb and a scaffold L50 of 7. We annotated a comprehensive set of genes that have been associated with models of aggressive behavior, a trait which the wolverine is purported to have in the popular literature. To support an integrated, genomics-based wildlife management strategy at a time of environmental disruption from climate change, we annotated the principal genes of the innate immune system to provide a resource to study the wolverine's susceptibility to new infectious and parasitic diseases. As a resource, we annotated genes involved in the modality of infection by the coronaviruses, an important class of viral pathogens of growing concern as shown by the recent spillover infections by severe acute respiratory syndrome coronavirus-2 to naïve wildlife. Tabulation of heterozygous single nucleotide variants in our specimen revealed a heterozygosity level of 0.065%, indicating a relatively diverse genetic pool that would serve as a baseline for the genomics-based conservation of the wolverine, a rare cold-adapted carnivore now under threat.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland.

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Identification of a new species of Ixodes Latreille, 1795 (Acari: Ixodidae), parasite of hog badgers (Carnivora: Mustelidae) in China.

Ixodes hunanensis n. sp. (Acari: Ixodidae), is identified based on the morphological characteristics and molecular biological analyses of males and females ex hog badger, Arctonyx collaris Cuvier (Carnivora: Mustelidae) from China. Adults of this new species are similar to those of other species of the subgenus Pholeoixodes Schulze, 1942, from which they can be distinguished by the shape of basis capituli, development of cornua, size of porose areas, shape, and size of spurs on coxae and phylogenetic analyses of the cox1 and 16S rRNA sequences.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae.

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes and searching for species-specific structural variants. Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observed species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This argues strongly for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications.

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and -0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNASer(AGC) and tRNALys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Early-life seasonal, weather and social effects on telomere length in a wild mammal.

Early-life environmental conditions can provide a source of individual variation in life-history strategies and senescence patterns. Conditions experienced in early life can be quantified by measuring telomere length, which can act as a biomarker of survival probability in some species. Here, we investigate whether seasonal changes, weather conditions and group size are associated with early-life and/or early-adulthood telomere length in a wild population of European badgers (Meles meles). We found substantial intra-annual changes in telomere length during the first 3 years of life, where within-individual effects showed shorter telomere lengths in the winter following the first spring and a trend for longer telomere lengths in the second spring compared to the first winter. In terms of weather conditions, cubs born in warmer, wetter springs with low rainfall variability had longer early-life (3-12 months old) telomeres. Additionally, cubs born in groups with more cubs had marginally longer early-life telomeres, providing no evidence of resource constraint from cub competition. We also found that the positive association between early-life telomere length and cub survival probability remained when social and weather variables were included. Finally, after sexual maturity, in early adulthood (i.e., 12-36 months) we found no significant association between same-sex adult group size and telomere length (i.e., no effect of intrasexual competition). Overall, we show that controlling for seasonal effects, which are linked to food availability, is important in telomere length analyses, and that variation in telomere length in badgers reflects early-life conditions and also predicts first year cub survival.

Comparative transcriptome analysis revealed omnivorous adaptation of the small intestine of Melinae.

As the main digestive organ, the small intestine plays a vital role in the digestion of animals. At present, most of the research on animal feeding habits focuses on carnivores and herbivores. However, the mechanism of feeding and digestion in omnivores remains unclear. This study aims to reveal the molecular basis of the omnivorous adaptive evolution of Melinae by comparing the transcriptome of the small intestines of Asian Badgers (Meles leucurus) and Northern Hog Badgers (Arctonyx albogularis). We obtained high-quality small intestinal transcriptome data from these two species. Key genes and signalling pathways were analysed through Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and other databases. Research has mainly found that orthologous genes related to six enzymes have undergone adaptive evolution. In addition, the study also found three digestion-related pathways (cGMP-PKG, cAMP, and Hippo). They are related to the digestion and absorption of nutrients, the secretion of intestinal fluids, and the transport of food through the small intestine, which may help omnivorous animals adapt to an omnivorous diet. Our study provides insight into the adaptation of Melinae to omnivores and affords a valuable transcriptome resource for future research.

Prevalence of Campylobacter spp. in Raccoon Dogs and Badgers in Miyazaki Prefecture, Japan.

A total of 55 samples of intestinal contents from 28 raccoon dogs (Nyctereutes procyonoides) and 27 badgers (Males anakuma) in Miyazaki prefecture, Japan, were examined for the presence of Campylobacter species. C. jejuni and C. upsaliensis were isolated from 3.6% (n = 1) and 75% (n = 21) of raccoon dogs, respectively. In contrast, no Campylobacter spp. was isolated from the badgers examined. The C. upsaliensis isolates were subjected to antimicrobial susceptibility testing against 8 antimicrobial agents. This revealed that most of the isolates from raccoon dogs were susceptible to the antimicrobial agents examined, whereas strains isolated from healthy dogs in Miyazaki prefecture, showed high rates of resistance. Virulence genes (flaA, cadF, ciaB, cdtA, cdtB, and cdtC) were present in the C. jejuni isolate from a raccoon dog, with the exception of flaB. By contrast, all these virulence genes examined were present in all C. upsaliensis strains isolated from raccoon dogs and dogs. The genetic diversity of those isolates based on the nucleotide sequences of 7 housekeeping genes (adk, aspA, atpA, glnA, glyA, pgi, tkt) was compared with that of C. upsaliensis strains isolated from dogs and strains selected randomly from humans and dogs deposited in the Campylobacter MLST database. The major cluster of raccoon dog strains was separated from both human and dog strains by phylogenetic tree analysis. These results suggest that raccoon dogs are a reservoir of C. upsaliensis and that isolates may represent a population different from that in humans and dogs. To our knowledge, this is the first study to have demonstrated a high prevalence of C. upsaliensis in raccoon dogs.

Genetic, social and maternal contributions to Mycobacterium bovis infection status in European badgers (Meles meles).

Within host populations, individuals can vary in their susceptibility to infections and in the severity and progression of disease once infected. Though mediated through differences in behaviour, resistance or tolerance, variation in disease outcomes ultimately stems from genetic and environmental (including social) factors. Despite obvious implications for the evolutionary, ecological and epidemiological dynamics of disease traits, the relative importance of these factors has rarely been quantified in naturally infected wild animal hosts. Here, we use a long-term capture-mark-recapture study of group-living European badgers (Meles meles) to characterize genetic and environmental sources of variation in host infection status by Mycobacterium bovis, the causative agent of bovine tuberculosis (bTB). We find that genetic factors contribute to M. bovis infection status, whether measured over a lifetime or across repeated captures. In the latter case, the heritability (h2 ) of infection status is close to zero in cubs and yearlings but increases in adulthood. Overall, environmental influences arising from a combination of social group membership (defined in time and space) and maternal effects appear to be more important than genetic factors. Thus, while genes do contribute to among-individual variation, they play a comparatively minor role, meaning that rapid evolution of host defences under parasite-mediated selection is unlikely (especially if selection is on young animals where h2 is lowest). Conversely, our results lend further support to the view that social and early-life environments are important drivers of the dynamics of bTB infection in badger populations specifically, and of disease traits in wild hosts more generally.

Genome Analysis of Sable Fur Color Links a Lightened Pigmentation Phenotype to a Frameshift Variant in the Tyrosinase-Related Protein 1 Gene.

Sable (Martes zibellina) is one of the most valuable species of fur animals. Wild-type sable fur color varies from sandy-yellow to black. Farm breeding and 90 years of directional selection have resulted in a generation of several sable breeds with a completely black coat color. In 2005, an unusually chocolate (pastel) puppy was born in the Puschkinsky State Fur Farm (Russia). We established that the pastel phenotype was inherited as a Mendelian autosomal recessive trait. We performed whole-genome sequencing of the sables with pastel fur color and identified a frameshift variant in the gene encoding membrane-bound tyrosinase-like enzyme (TYRP1). TYRP1 is involved in the stability of the tyrosinase enzyme and participates in the synthesis of eumelanin. These data represent the first reported variant linked to fur color in sables and reveal the molecular genetic basis for pastel color pigmentation. These data are also useful for tracking economically valuable fur traits in sable breeding programs.

An Evolutionary Cancer Epigenetic Approach Revealed DNA Hypermethylation of Ultra-Conserved Non-Coding Elements in Squamous Cell Carcinoma of Different Mammalian Species.

BACKGROUND: Ultra-conserved non-coding elements (UCNEs) are genomic sequences that exhibit > 95% sequence identity between humans, mammals, birds, reptiles, and fish. Recent findings reported their functional role in cancer. The aim of this study was to evaluate the DNA methylation modifications of UNCEs in squamous cell carcinoma (SCC) from different mammal species. METHODS: Fifty SCCs from 26 humans, 17 cats, 3 dogs, 1 horse, 1 bovine, 1 badger, and 1 porcupine were investigated. Fourteen feline stomatitis and normal samples from 36 healthy human donors, 7 cats, 5 dogs, 5 horses, 2 bovines and 1 badger were collected as normal controls. Bisulfite next generation sequencing evaluated the DNA methylation level from seven UCNEs (uc.160, uc.283, uc.416, uc.339, uc.270, uc.299, and uc.328). RESULTS: 57/59 CpGs were significantly different according to the Kruskal-Wallis test (p < 0.05) comparing normal samples with SCC. A common DNA hypermethylation pattern was observed in SCCs from all the species evaluated in this study, with an increasing trend of hypermethylation starting from normal mucosa, through stomatitis to SCC. CONCLUSIONS: Our findings indicate that UCNEs are hypermethylated in human SCC, and this behavior is also conserved among different species of mammals.

Comparative Analysis of the Umami Taste Receptor Gene Tas1r1 in Mustelidae.

Pseudogenization of the umami taste receptor gene Tas1r1 has been reported more commonly in aquatic than terrestrial mammals. We hypothesized that the more species are adapted to the aquatic environment, the less important a role the umami taste detection plays. To test this hypothesis, we focused on the Mustelidae family because their habitat and feeding ecology are highly diverse. We found pseudogenizing mutations in exon 1 of Eurasian otter and exon 6 of African clawless otter, both of which forage underwater. Our analysis of the ratio of nonsynonymous to synonymous nucleotide substitution rates suggested that purifying selection pressures on Tas1r1 are weaker in the lineages with non-functional Tas1r1 than the lineages retaining functional Tas1r1. Our analysis also suggested that relaxed selection pressures on Tas1r1 in Mustelidae species adapted to the aquatic environment, although we cannot exclude the possibility that they are restricted to Lutrinae irrespective of their feeding habitat. Overall, the results of the present study support the idea that differences in selection pressures on Tas1r1 reflect differences in feeding behaviors.